Finite-size scaling for the Ising model on the Mobius strip and the Klein bottle 
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We study the finite-size scaling properties of the Ising model on the Mobius strip and the Klein 
bottle. The results are compared with those of the Ising model under different boundary conditions, 
that is, the free, cylindrical, and toroidal boundary conditions. The difference in the magnetization 
distribution function p(m) for various boundary conditions is discussed in terms of the number of 
the percolating clusters and the cluster size. We also find interesting aspect-ratio dependence of the 
value of the Binder parameter at T — T c for various boundary conditions. We discuss the relation 
to the finite-size correction calculations for the dimer statistics. 



PACS numbers: 75.10.Hk, 05.50.+q, 64.60.-i 

The systems under various boundary conditions have 
the same per-site free energy in the bulk limit, whereas 
the finite-size corrections are different. The idea of finite- 
size scaling (FSS) p],^] is very important in understand- 
ing finite-size effects near the criticality. It is known that 
the FSS functions depend on the boundary conditions. 
The difference in the FSS functions for the Ising model 
under the periodic and free boundary conditions has been 
discussed in connection with the universal FSS for the 
percolation problem || and the Ising model Q . The sys- 
tems with the tilted boundary conditions have been also 
studied J|||. Quite recently, Lu and Wu have stud- 
ied dimer statistics on the Mobius strip and the Klein 
bottle. These two systems are other examples of inter- 
esting boundary conditions, and their topological prop- 
erty is unique. The dimer statistics on the Mobius strip 
has been also studied by Brankov and Priezzhev || . The 
ground-state entropy of the Potts antiferromagnet on the 
Mobius strip has been investigated |j] . In two dimensions 
(2D) the relevance of the finite-size properties to the con- 
formal field theory is another source of interest |K]-n2] . 
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FIG. 1. Illustration of the rectangular lattice with various 
boundary conditions; plane, cylinder, torus, Mobius strip, and 
Klein bottle. Thick and thin lines denote periodic and free 
boundaries, respectively. 



In this Letter we study the FSS functions for the 2D 
Ising model on the Mobius strip and the Klein bottle in 
view of increasing interest on the effect of boundary con- 
ditions for finite systems. We compare these results with 



those of the Ising model under different boundary condi- 
tions, that is, the free, cylindrical, and toroidal boundary 
conditions. As a total five boundary conditions are con- 
sidered, and they are illustrated in Fig. |l|. We deal with 
the rectangular lattice of size Li x L 2 with the aspect 
ratio a = L\jLi- Thick and thin lines denote periodic 
and free boundaries, respectively, in Fig. [l| We impose 
the periodic boundary conditions in both the horizontal 
and vertical directions for the torus (toroidal boundary 
condition). The twisted periodic boundary condition is 
imposed in the y direction, that is, 



f(x,y + L 2 ) = f{Li ~x,y), 



(1) 



for the Klein bottle. The periodic boundary condition is 
imposed in one direction and the free one in the other di- 
rection for the cylinder (cylindrical boundary condition); 
the twisted periodic and the free boundary conditions 
are imposed for the Mobius strip. We refer to our sys- 
tem as a plane when we impose the free boundary condi- 
tions in both directions. The Mobius strip and the Klein 
bottle have a unique topological property; they have a 
non-orientable surface. The topological properties of five 
boundary conditions, the number of sides and the num- 
ber of edges, are tabulated in Table |. The symmetry 
property under a transformation a — > 1/a is also given in 
Table [I]. The system is symmetric under such a transfor- 
mation for the plane and the torus. 



TABLE I. The topological properties of the rectangular 
lattices with various geometries (boundary conditions); the 
number of sides and the number of edge are given. The sym- 
metry property under the transformation a — > 1/a is also 
given. 



geometry 



# of side 



# of edge 



1/a 



Klein bottle 
Mobius strip 

torus 
cylinder 

plane 



no 
no 

yes 
no 
yes 
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We use the Monte Carlo simulation to study the FSS 
properties of the 2D Ising model with various boundary 
conditions near the criticality. The moments of magneti- 
zation are basic quantities for the FSS analysis. Here we 
focus on the Binder parameter 13 1, which is defined by 



(m 4 ) 



Vs- 

2 V (m 2 ) 2 



(2) 



One may determine the critical point from the crossing 
point of g for different sizes as far as the corrections to 
FSS are negligible. The value of the Binder parameter at 
the critical point is not a universal quantity, and it de- 
pends on the shape of the finite systems and the bound- 
ary conditions. 



0.5 



closed marks 80x20 
open marks 160x40 




1 

mL 0125 



0.5 



K 



*«^ T '•a ^ *s.. 



closed marks 80x20 
open marks 1 60x40 



-3 



-2 



-1 



1 

(T-TJ L 



FIG. 2. Plot of g as a function of (T-T c )L 1/u for the Ising 
model with various boundary conditions, where T c — 2.269 ■ ■ ■ 
and u=l. The data are plotted in units of the coupling J. 
Here, K, M, T, C, and P represent Klein bottle, Mobius strip, 
torus, cylinder, and plane, respectively. The aspect ratio a is 
chosen as 4. 



We plot the temperature dependence of the Binder pa- 
rameter g for several lattices with different boundary con- 
ditions in Fig. |[ The system sizes are given within the 
figure, and the aspect ratio of the lattice is chosen as 
a = 4. The data of different sizes collapse on a single 
curve with using the scaling variable (T — T C )L 1 ^ V for the 
horizontal axis. Here, L = (L x x L 2 ) 1/2 , T c = 2.269- • ■ 
and v = 1 for the 2D Ising model; we plot the data in 
units of the coupling J. We have very good FSS behavior 
for g, and also find the strong dependence on the bound- 
ary conditions. The g values of the Klein bottle are larger 
than those of the torus. The same behavior is found in 
the g values of the Mobius strip and those of the cylinder. 



FIG. 3. Plot of p(m)L _,3/ ' / at T = T c as a function of 
m jJ } l v f or the Ising model with various boundary conditions, 
where /3/u=l/8. Here, K, M, T, C, and P represent Klein 
bottle, Mobius strip, torus, cylinder, and plane, respectively. 
The aspect ratio a is chosen as 4. 

In order to clarify the difference in the Binder param- 
eter g, we study the magnetization distribution function 
p{m) at T = T c . We show the FSS plots of p(m) for vari- 
ous sizes with different boundary conditions in Fig. |^. We 
plot p(m)L- /3/,y as a function of mL^l v \ (3 = 1/8 for the 
2D Ising model. We again have good FSS behavior for 
p(m) at T = T c . The FSS function of pirn) strongly de- 
pends on the boundary conditions. There are two sharp 
peaks in p(m) for the Klein bottle and Mobius strip, and 
the secondary broad peak around m = becomes larger 
for the torus and cylinder. There is one peak around 
to = for the plane. 

We may understand this behavior from the number of 
percolating clusters and the size of clusters |14|. For the 
systems with large aspect ratio, the importance of the 
number of percolating clusters has been pointed out by 
Hu and Lin |]l5f . The probability for the appearance of n 
percolating clusters for anisotropic lattices have been of 
current interest p6|^7[ . It has been revealed by Tomita 
et al [14[ that for the Ising model the combination of the 
percolating clusters with up spins and those with down 
spins gives the contribution to the broad peak around 
777 = in p{m). It is interesting to relate the effect of 
the boundary conditions to the number of percolating 
clusters and the size of clusters. The periodic bound- 
ary conditions may have the tendency that the size of 
the percolating clusters becomes larger compared to the 
free boundary conditions; the order will be reduced near 
a free surface. Moreover, if one twists the boundaries, 
more clusters mix together. To confirm this speculation, 
we study the percolating properties of the Ising model 
with different boundary conditions. Using the fact that 
the Ising model is mapped to the percolation problem 
with the bond concentration of 1 — e~ 2J / T [[18,191, we 
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can assign clusters. Then, we may decompose the physi- 
cal quantities by the number of percolating clusters [fL4| . 
In Fig. H we plot the fraction of lattice sites in the n per- 
colating clusters (c}„ at T = T c for the system of size 
160 x 40. We find from the figure that the cluster sizes 
become smaller for the system with free bondary condi- 
tions, and the single cluster is actually more dominant 
for the Klein bottle and the Mobius strip. Thus, we ex- 
plain the large g values of the Mobius strip and the Klein 
bottle in relation to the number of percolating clusters. 
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for a < 1. For small enough a (<C 1), the Mobius strip, 
the cylinder, and the plane show the same FSS proper- 
ties because the boundaries along the shorter directions 
for these three are the same, that is, the free boundary 
condition. 
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FIG. 5. Aspect ratio dependence of the Binder parameter 
at T — T c , g c , for the Ising model with various boundary con- 
ditions; that is, Klein bottle, Mobius strip, torus, cylinder, 
and plane. 



FIG. 4. Plot of (c) n at T — T c as a function of n for the 
Ising model with various boundary conditions; that is, Klein 
bottle, Mobius strip, torus, cylinder, and plane. The system 
size is 160 x 40 (a = 4). 

The FSS functions for the Binder parameter g and the 
distribution function p(m) depend on both the aspect ra- 
tio and the boundary condition. Wc plot the aspect-ratio 
dependence of g at T = T c , g c , for various boundary con- 
ditions in Fig. H We use the logarithmic scales for the 
horizontal axis. We have several interesting observations 
from Fig. |^. The plot of g c as a function of a in log- 
arithmic scale is symmetric for the torus and the plane 
because of the symmetric property under the transforma- 
tion a — > 1 /a; g c takes the maximum at a = 1. This is not 
the case for other geometries. For large enough a (>■ 1), 
the FSS properties of the Klein bottle and those of the 
Mobius strip become the same because the boundaries 
along the shorter direction determine the FSS properties 
of the system; for both the Mobius strip and the Klein 
bottle, the boundary condition along y axis is the twisted 
periodic one. The FSS properties of the torus and the 
cylinder are the same for large enough a. In contrast, 
the systems with a < 1, the Klein bottle and the torus 
show similar FSS behavior. It is because the twisted pe- 
riodic boundary or simple periodic one is not important 
for the number of percolating clusters for a < 1. It is 
the same situation for the Mobius strip and the cylinder 



Let us compare our results with those of the dimer 
statistics by Lu and Wu 0, who have calculated the 
finite-size corrections for the dimer generating function. 
The finite-size correction coefficients have been tabulated 
in Table 1 of Ref. @. Their notation M x N for the sys- 
tem size corresponds to our Li x ^2- In the limit of 
large M in their notation, which corresponds to a > 1, 
the finite-size correction coefficients (C2, A2) are classified 
into three groups, that is, Mobius strip and Klein bottle, 
cylinder and torus, and plane. In contrast, for large N 
in their notation, a <C 1, the finite-size correction coeffi- 
cients (ci,Ai) are classified into two groups; one group 
is Klein bottle and torus, and the other is Mobius strip, 
cylinder, and plane. For both cases their results of the 
dimer statistics are consistent with the present results. In 
addition, for large M = N, which corresponds to a = 1, 
the correction coefficients (ci, C2) are relevant; then Klein 
bottle and torus, Mobius strip and cylinder, and plane 
form three groups. This result is again consistent with 
the present one. These observations are compatible with 
the conformal field theory [|l~0[-p~2[| that the behavior of the 
difference in finite-size corrections for different boundary 
conditions are model- independent. 

To summarize, we have studied the FSS properties of 
the Ising model with various boundary conditions, that 
is, Klein bottle, Mobius strip, torus, cylinder, and plane. 
We have elucidated the difference in the magnetization 
distribution function p(m) for various boundary condi- 
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tions in terms of the number of the percolating clusters 
and the cluster size. We have found interesting aspect- 
ratio dependence of the g value at T — T c for various 
boundary conditions. The FSS properties of the systems 
are classified into three groups for s>l and two groups 
for a <C 1. For large a, the g value becomes large for the 
Mobius strip and the Klein bottle, which is characteristic 
for systems with a non-orientable surface. 

There may be several directions for future study. In 
the dimer statistics, two identities relating dimer gen- 
erating functions for Mobius strips and cylinders have 
been established 0. It is desirable to explore exact re- 
lations between the Binder parameters of the Ising sys- 
tems with different boundary conditions. Anyons on the 
cylinder and the torus have been studied with the braid- 
group analysis p0| ; the topology of the systems is im- 
portant in the fractional quantum Hall effect. Anyons on 
the Mobius strip and the Klein bottle will be interesting 
subjects to study. The boundary condition dependence 
of the critical behavior of the Anderson transition has 
been investigated (22) . It is also interesting to study the 
scaling functions for the Anderson transition with the 
Mobius and Klein-bottle boundary conditions. 

Several systems on the Mobius strip and the Klein bot- 
tle have been studied in various fields of physics. Persis- 
tent currents in a Mobius ladder have been studied [^l| , 
and interesting finite-size effects have been discussed. A 
general construction of correlation functions in rational 
conformal field theory on the Mobius strip and the Klein 
bottle has been made in terms of three-dimensional topo- 
logical quantum field theory pq |. Moreover, the matrix 
string theory has been constructed on the Mobius strip 
and the Klein bottle |24j] . The present work on the inter- 
play of the topology and aspect ratio in the FSS proper- 
ties of the Ising model may accelerate the study of physics 
on a non-orientable surface. 

After we submitted our manuscript, we came to know 
the preprint by Lu and Wu |2^| . They have studied the 
partition function of the Ising model on the Mobius strip 
and the Klein bottle analytically. They obtained essen- 
tially the same conclusion as that of the dimer statistics 
. We have made the detailed study of the boundary- 
condition dependence of the magnetization distribution 
function p(m) , which is complementary to the analytical 
work on the partition function. 
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